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ABSTRACT 

o, 

rj ^ We present the global results of a large spectroscopic survey carried out in order to 

identify z ~ 5 Lyman break galaxies (LBGs) across ten widely-separated ^ 45 arcmin'^ 
fields to a depth oi Iab = 26.3. The redshifts of seventy 4.6 < z < 5.6 LBGs were iden- 
tified through their Lya emission and/or a strong continuum break, with thirty eight 
Q ■ sources showing detectable line emission of between 2.6 x f 0~^^ and 7 x 10~^^ erg cm~^ 

H ' s~^. Just over half of the spectroscopically-confirmed z ^ 5 galaxies have rest-frame 

C/3 , Lya equivalent widths above 20 A, double the frequency of similarly strong line emit- 

^, ' ters in similar z ^ 3 LEG samples. However, when reasonable corrections are made for 

the spectroscopically-unconfirmed sources that are nevertheless at these redshifts in 
both samples, we find no significant difference in the frequency of high equivalent- width 
^ ' line emitters between the samples. The rest-frame UV continuum slope of a typical 

t~^ I z ~ 5 line-emitting galaxy (as measured primarily from photometry, but also apparent 

^^ I in spectroscopy) is bluer than that of a typical break-only galaxy, a difference that is 

OO ' difficult to explain purely by differences in the ages of their stellar populations. Varia- 

^^ , tion in metallicity and/or dust extinction can more straightforwardly account for this 

[■"-« ■ difference. If the correlation between metallicity and UV continuum slope identified at 

^^ ' low redshift is applicable at z > 3, the typical z ^ 5 LBGs have metallicities a factor 

^D . of three lower than those of LBGs at z ^^ 3. HST imaging of a subset of the LBGs 

indicates that a large majority of the spectroscopically-confirmed LBGs in our sample 
are members of multiple systems (on ~ arcsec scales) and/or show disturbed morphol- 
ogy. Using local LBG analogues as a model, this multiplicity could be explained either 
^^ ' by super-starburst regions within a larger unseen structure, or by a high incidence of 

JH I merging events at this epoch. The current data cannot distinguish between these two 

possibilities. The surface density of z ~ 5 LBGs in two of the ten fields is considerably 
higher than in the rest. Both show clear spikes in their redshift distributions indicating 
strong three-dimensional clustering in these fields. Against an expectation of about 
one source per 0.1 in redshift between 4.8 < z < 5.6, one field has seven identified 
objects between 5.11 < z < 5.21 and the other has 17 between 4.95 < z < 5.15. Nei- 
ther structure can be bound given their depth in redshift and probably extend beyond 
the observed fields. The three-dimensional distances between LBGs in the structures 
are too large for them to have triggered their starbursts through mutual gravitational 
interaction, and so it is likely that the short-lived LBGs represent only a small fraction 
of the baryons in the structures. 
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1 INTRODUCTION plex physical processes can affect the properties of galaxies. 

Consequently, it is hoped that the earlier a galaxy is ob- 



The identification and observational study of the earliest 
galaxies is crucial to our understanding of how galaxies form 
and evolve. With increasing cosmic age, increasingly com- 



served, the more straightforward it will be to compare its 
properties to any theory of galaxy growth. Whether or not 
this turns out to be true in practice, such galaxies can be 



2 L. S. Douglas et al. 



used as markers of early structure formation and explored 
as candidate emitters of the radiation responsible for reion- 
ization. 

While several different techniques have been success- 
fully used to identify galaxies at z > 5, the Lyman break 
technique is arguably the most useful as it can provide 
a sample of objects with measurable properties for indi- 
vidual galaxies (Lyman a line strength, continuum spec- 
tral energy distribution and potentially absorption line di- 
agnostics of the sample or of individual objects if bright 
enough). Optically-selected samples of high redshift Lyman 
break galaxies (LBGs) are dominated by systems that are 
undergoing strong and unobscured star formation because 
by 2 ~ 5 optical observations probe the rest-frame spectral 
region shortward of I8OOA, where emission is dominated by 
comparatively short-lived O and B stars. 

In this decade there has been considerable progress in 
the observational identification of systems at these redshifts. 
This has allowed us to take the first steps towards under- 
standing the properties of galaxies that most likely represent 
the earliest st ages in the evoluti on of current-era massive 
galaxies ( e.g. IVerma et al.l 120071 ). Thousands of candidate 
2; ~ 5 star forming LBGs have been identified in large- 
area ground-based photometric surveys, but only a very 
small fraction of these have been confirmed spectroscopi- 
cally. While studies of optically-selected photometric sam- 
ples can provide reliable insights into the properties of the 
sources when HST imaging is used on its o wn or in combi- 
nation with deep infra-r ed photometry fe.q. iBouwens et al.l 
I2OO7I : IVerma et al.ll2007f ). deriving reliable conclusions from 
ground-based optically-selected samples is more problem- 
atic. Not only do thes e samples suffer from the effects of 
conta mination (see e.g. IStanwav et al.ll2008bl : lDouglas et al.l 
I2OO9I . hereafter paper I), but the information supplied on 
each system is minimal and so their detailed properties can- 
not be explored. 

Spectroscopically-confirmed samples of z ~ 5, while in- 
evitably smaller than photometric samples, ultimately allow 
more detailed studies of the properties of the galaxies. With 
a reliable sample, uncertainties in statistical properties such 
as the luminosity function can be minimised. Correlations 
between observables, such as line-strength and colour, start 
to unravel the evolution of the objects. Instead of explor- 
ing the clustering by statistical methods in two dimensions 
(which can be unreliable if the sample is contaminated), a 
direct determination of three-dimensional clustering can be 
made. With the stochastic and relatively short-lived nature 
of the star formation within these systems, the interpreta- 
tion of clustering is potentially complicated given that the 
timescale over which a source may be UV-luminous can be 
only a small fraction of the look-back time through a survey 
volume. The most important benefit of a spectroscopically- 
confirmed sample is that it allows further spectroscopi c 
follow-up in other wavebands {e.g. IStanwav et al.ll2008al ). 
so will eventually lead to a fuller understanding of the early 
stages of galaxy formation and evolution. 

In order to furnish ourselves with a sample of objects 
with which to carry out such studies, we have carried out a 
spectroscopic survey of ten widely-separated ~ 45 arcmin^ 
fields as part of the ESO Remote Galaxy Survey (ERGS), 
an ESO large programme (PI M. Bremer, ID 175.A-0706) in 
order to identify and spectroscopically-confirm a sample of 



z ~ 5 unobscured star-forming galaxies. We had several mo- 
tivations for carrying out this spectroscopic survey. Firstly, 
to explore the reliability of photometric samples and how 
they impact the reliability of statistical properties for 2 ~ 5 
LBGs, e.g. their luminosity function (see paper I). Secondly, 
to obtain as large a sample of confirmed z ~ 5 LBGs as pos- 
sible in a reasonable amount of time (the aim being at least 
50 2: ~ 5 LBGs in 100 hours of 8m spectroscopy). A sample 
of this size (especially when combined with other published 
samples) allows exploration of the properties of the sources 
and of any correlation between them. Thirdly, by obtain- 
ing such a sample over multiple widely-separated pointings, 
3-dimensional clustering can be explored directly, without 
being dominated by cosmic variance in any one field. 

The structure of the paper is as follows. In section [5] 
we describe the data and sample selection of high redshift 
LBGs. Details and examples of spectroscopic identification 
are given in section [31 Discussion of the global properties of 
the sample follows in subsequent sections. The z ~ 5 redshift 
distribution is described in section U Section [5] explores the 
rest-frame UV continuum of the LBGs and in section [S] the 
properties of the observed Lya emission are discussed. In 
section [7| we describe the morphology of the confirmed high 
redshift LBGs. The surface density of sources is discussed 
in section [8] and the results of the z ~ 6 LBG search are 
presented in section [O] 

We adopt the standard ACDM cosmology, i.e. a flat uni- 
verse with Oa = 0.7, Q.M = 0.3 and Ho = 70km s"^ Mpc"^ 
All p hotometry was det ermined in the AB magnitude sys- 
tem (|Oke fc Gunnll 19831 ). 



2 DATA AND SAMPLE SELECTION 

2.1 Imaging Data 

The spectroscopic sample was derived from the same set 
of imaging data used to derive the photometric sample 
previously presented in paper I. The data included deep 
(around two hours per band) V , R, I and z-b and imaging 
obtained with FORS2 (|Appenzeller et al.ll 199^ 1 on the VLT 
and moderately-deep J and Ks-hand imaging (5-6 hours per 
band) obtained with SOFI on the NTT. The optical data 
covered approximately 49 arcmin^ in each field, while the 
near-IR data covered approximately half the area of the 
optical data per field. While the photometric sample re- 
ported in paper I was drawn from 286 arcmin^ requiring 
both optical and near-IR coverage, the spectroscopic candi- 
dates were drawn from an effective area of 440 arcmin^ , the 
area covered by the R and 7-band imaging after taking into 
account the area compromised by high surface brightness 
foreground objects. All imaging except the z-band frame 
were obtained as part of the ESO Distant Gluster Survey 
(EDisCS), itself a n ESO large prog ramme (PI S. White, 
ID:166.A-0162, see I White et al.ll2005l ). The ZAB-h&nd imag- 
ing was obtained as part of ERGS. Although there was 
some field-to-field variation, typical 2a limits in a 2" circular 
aperture were yAs=28.1, Rab=28.0, Iab=27.1, 2as=26.0, 
Jas=24.5, Kab=23.6 and seeing varied between ~ 0.6" and 
0.8" in the optical bands, 0.8" and 1.0" in Jab and 0.6" and 
0.8" in Kab, comparable to the quality and depth of other 
surveys of the equivalent or larger area. See paper I for a 
detailed field-to-field breakdown. 
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2.2 Spectroscopic Sample 

As for the photometric sample discussed in paper I, the pri- 
mary selection criteria of Iab < 26.3, Rab — Iab > 1-3 
and Vab < 27.5 were used to select objects that may have 
the characte ristic spectral break of z > 4.8 LBGs. The same 
SExtractor (jBertin fc ArnoutdllQQg ') catalogues were used 
to select objects for the photometric sample and spectro- 
scopic targets. The photometric sample was generated us- 
ing 2 arcsec diameter circular apertures so that photometry 
across all the bands (with different pixel sizes) could be com- 
bined in the object selection. While the photometric sample 
was defined to allow an accurate determination of complete- 
ness and contamination, one aim of the spectroscopy was 
to maximise the number of spectroscopic candidates to be 
targeted with the slit masks, thereby maximising the size 
of the spectroscopically-confirmed z ^ 5 LBG sample. Con- 
sequently, we performed the candidate selection twice using 
the same criteria with different photometric catalogues. The 
first catalogue used the aperture magnitudes from paper I 
and the second used total magnitudes measured using the 
SExtractor MAG_AUTO apertures defined in I and with ap- 
propriate corrections to the V and _R-band value (typically 
of a few hundredths of a magnitude due to the difference in 
seeing between the bands). If an object satisfied the primary 
selection criteria in either of the apertures, it was flagged as 
a candidate for spectroscopy. In addition, as we specifically 
wanted to both maximise the number of candidates placed 
on the spectroscopic masks and explore the potentially con- 
taminating populations, we avoided an explicit Iab — zab 
colour cut to rule out potential lower redshift interlopers. In 
designing the observational strategy, we matched the num- 
ber of masks to be used in any field to the number of ob- 
jects selected by the primary selection criteria, ensuring that 
there were generally sufficient slits to observe objects with 
any Iab — zab colour. Figure [T] shows the range of colours 
of the objects observed during the spectroscopy program. 

Even with several masks per field, pairs of candidates 
will compete for space because their positions are such that 
their dispersed spectra would overlap. Consequently the can- 
didates were prioritised drawing on our previous experience 
of comparable spectroscopy with the same instrumentation. 
We had previously observ ed four adjoining FORS 2 fields, 
searching for z ~ 5 LBGs (jLehnert fc Bremen 120031 . Davies 
et al. in prep). As a result of this experience we were able to 
prioritise the candidates based on their likelihood of being 
at high redshift, as described below. This was done to max- 
imise the number confirmed of 2 ~ 5 LBGs on the masks. 

Each candidate was categorised according to its pho- 
tometry in all bands and if necessary from visual inspec- 
tion of all of the ground-based images. If the object had 
Iab — Zab < 1, was undetected in the V, J and _R's-band 
images, and the photometry was in no way influenced by 
neighbouring objects (affecting the background estimate or 
the positioning of the aperture), the object was believed to 
be at high redshift and given the highest priority (noted as 
priority 1 in Tables 1 and 2). If the object was detected in the 
V-band, indicating that it may be a lower redshift galaxy or 
a star, or had potentially compromised photometry due to 
object overcrowding it was assigned a lower priority (prior- 
ity 2 in Tables 1 and 2). Sources which were detected in the 
near-infrared imaging, with colours inconsistent with high 
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Figure 1. Colour-colour diagram of the spcctroscopically ob- 
served objects illustrating the range of colours sampled. The 
dashed line is the formal colour cut used in the photometric se- 
lection. Arrows depict where 2a limits were used. 



redshift star-forming galaxies were classified as likely 'low z' 
objects. Finally, in empty areas of the spectroscopic masks 
where no other strong candidates could be placed, objects 
with photometry outside the selection criteria {e.g. objects 
with red Rab — Iab colours but Iab >26.3, objects with 
1.0 < Rab — Iab < 1.3, or objects detected in the l/-band) 
were observed, classified as 'other' sources. This was partic- 
ularly valuable as it probed additional colour space outside 
the selection. In practice, sufficient objects were observed to 
obtain good statistics for each class despite slit contention. 

Our previous experience had shown that brighter 
{Iab < 25) candidates were almost invariably low redshift 
galaxies (0.6 < z < 1.2) or stars, especially if they had red 
Iab — Zab or Iab — Kab colours (see paper I for details 
of the colour modelling of the contaminating populations) . 
However, even in the cases where the multi-band character- 
istics of an object strongly indicated a low or zero-redshift, 
an object was not rejected from the list of candidates, merely 
given a lower priority. This was done so that the reliability 
of our selection, both for this sample and also for the photo- 
metric sample in paper I, could be determined. For objects 
with Iab < 25, it is impossible to rule out with 100 per 
cent certainty an object being at high redshift based on the 
ground-based photometry alone; all that can be determined 
is that an object is highly unlikely to be at high redshift. 
Objects brighter than Iab < 24 selected by our colour cuts 
would most likely have a strong AGN component in order 
to be at 2 ~ 5. While our optical colours do not distinguish 
between a z ~ 5 quasar and certain types of late-type stars, 
given that we detected only a single fainter object with a 
mode rate AGN component in our sample (see lDouglas et al.l 
[20071) , it is likely that no brighter AGN existed in the volume 
probed by the survey. 

As the near-IR data did not cover all of the optical 
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fields, the full prioritisation could not be carried out uni- 
formly across all areas. This resulted in less rigorous pri- 
oritisation in a fraction of each field. Fortunately, objects 
which were near-IR detected and therefore less likely to be 
at the highest redshifts tended to be the brighter objects in 
the optical, or had red Iab — zab colours, both of which fac- 
tors downgraded their priority on their own. Thus a lack of 
near-IR data did not dramatically change the prioritisation 
of objects between the regions with and without near-IR 
data. As prioritisation was only important when there was 
contention between candidates for space on the slit masks, 
many candidates were observed regardless of their priority. 
Table [T] shows the fraction of high and low redshift targets 
which were observed. The results in table[2] (discussed below) 
confirm a posteriori that the prioritisation was effective. 

In addition to the 2 ~ 5 target selection, our data 
allowed us to select z ~ 6 candidates using a redder 
Iab — ZAB > 1.3 col our cut which at zab > 24.5 will select 
LBGs at z > 6 (e.g. IStanwav et aLlbOO^ ) . Such objects will 
not be detected in our near-IR imaging, or in the V-band im- 
age. These higher redshift LBG candidates were rarer than 
those at z ~ 5. There were typically nine such candidates 
per field and a total of 73 candidates were observed on the 
masks, placed in regions where they did not compete for 
space with the highest priority 2 ---^ 5 candidates. 

2.3 Observations 

Twenty VLT/FORS2 'MXU' spectroscopic masks were used 
to obtain redshifts in the ten fields. The number of masks 
assigned to each field broadly refiected the number of high 
priority z -^ 5 candidates. The only exception is that the 
field of J1103.7 received an extra mask at the expense of 
field J1216.8 because the placement of candidates in the 
latter field made for an inefficient third mask. Each mask 
was designed to have between 30 to 40 slitlets of width one 
arcsecond and typical length of 10 arcsec. In general, where 
two or more candidates contended for space on a slit mask, 
the highest priority object was chosen. Occasionally, in order 
to help calibrate the reliability of the photometric sample, a 
potential contaminant was chosen at random over a higher 
priority object. For this reason, and because high priority 
candidates will inevitably contend with each other for space 
on the mask, not all of the high priority targets had slitlets 
assigned to them (see table [TJ. 

The observations with FORS2 were carried out in ser- 
vice mode on nights between December 2005 and May 2006. 
Seeing was always better than 1 arcsec and observations were 
taken during dark time, clear or photometric conditions and 
at an airmass below two. The masks were made using fims, 
the mask design software for F0RS2 distributed by ESO. 
Each mask was observed in five blocks of 43 minutes. Each 
block consisted of four observations of 650 seconds at four 
different nod positions along the slit resulting in a total on 
source exposure time of 3.6 hours. These nod positions were 
also offset between blocks ensuring that the spectrum of 
each object was placed at 20 independent positions on the 
CCD. This technique was used to avoid the effects of pat- 
tern noise or bad pixels and to facilitate the best possible 
sky subtraction. FORS2 was used in spectroscopic mode us- 
ing the 3001 grism with the OG525 blocking filter to restrict 
the wavelength range to between 650-lOOOnm. The resulting 



pixel scale was 3.2A in the spectral direction by 0.25 arcsec 
spatially and the spectral resolution was ~ lOA, measured 
from the width of the s ky lines. This observing s etup was 
identical to that used bv lLehnert fc Bremeij (|2003l ) and the 
choice of grism and b locki ng filter the same as t hat used by 
IVanzeUa etHI (120051 ) and lVanzella et al.l (|2009l ). 



2.4 Data Reduction 

The data were reduced foll owing the same standard proce- 
dures as those described in lLehnert fc Bremen (|2003h . with 
the following modifications. The individual spectra were rec- 
tified before sky subtraction. As the seeing of the spectro- 
scopic observations was close to (but always less than) the 
slit width, the spectroscopy did not recover all of the fiux 
of an object. The standard star spectra were taken through 
a 5 arcsec-wide slit and so recovered a higher fraction of 
the stellar fiux. To account for this, the ratio of flux re- 
covered in 1" and 5" wide rectangular apertures was mea- 
sured from Jyis-band imaging of standard stars, and the 
fluxes determined in the spectroscopy were corrected by 
the factor appropriate to the seeing of the individual ex- 
posures. The correction factors typically ranged from 1.1 
to 1.5. The final combined exposures reached a fiux limit 
of a few 10~^* ergscm"^ s~^ A~^ for an unresolved emis- 
sion line in a dark region of the sky spectrum, compara- 
ble to the depth obtained by iLehnert fc Bremen (|2003l ) as 
expected given the identical observational setup (exposure 
time, grism, slit width). 



3 SPECTROSCOPIC IDENTIFICATIONS 

The reduced spectra were examined by four of the authors 
(LSD, MNB, ERS and MDL) and classified into high red- 
shift, low or intermediate redshift and low signal-to-noise 
categories. The high redshift sources were expected to show 
a clear continuum break at the redshifted wavelength of 
Lya, often with a Lya emission line at the same wave- 
length. Low or intermediate redshift objects were identified 
through multiple emission lines (such as II/3 and [GUI]), no 
break in the continuum or a steeply reddening continuum 
with molecular absorp tion bands indicative of low mass stars 
(jStanwav et al.ll2008d ). Low redshift objects were also iden- 
tified by the presence of continuum emission at the bluest 
wavelengths (~6500A), inconsistent with a red Rab — Iab 
colour and therefore included in the sample due to photo- 
metric uncertainty. Any spectra exhibiting one or more of 
these features were classified as low redshift. The object did 
not necessarily have to have a secure redshift, merely that 
the spectrum was clearly inconsistent with that of a high 
redshift LBG. Objects where the signal-to-noise in the con- 
tinuum was too low to be placed in either of the above cat- 
egories were classified as 'unknown'. 

In addition to categorising an object as being at high 
redshift (z ~ 5), we split this category into two grades: A 
and B. If the spectrum was clearly that of a high redshift 
object with an emission line at the wavelength of a con- 
tinuum break, or a clear continuum break with no other 
obvious low redshift features red-ward or blue-ward of the 
break, the spectrum was classified as grade A. The redshift 
was determined to the resolution limit of the spectrograph. 
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Table 1. The number of targets which were observed spectroscopically in the full optical field of view. Targets include high redshift 
candidates divided into priority 1 and 2 (see text) and likely lower redshift candidates which obey the optical selection criteria but are 
detected in the near-IR imaging. For each field, the number of photometrically selected targets and the number observed on mask are 
given for each target category. 



Field 


No. 


Priority 


1 2 


:-5 


Priority 


2 z~5 


Low z 


tarj 


2;ets 




masks 


targets 


on 


mask 


targets 


on 


mask 


targets 


on 


mask 


J1037.9-f243 


f 


7 




5 


8 




2 


5 




3 


J1040.7-lf55 


3 


30 




24 


12 




4 


12 




7 


Jf054.4-lf46 


f 


7 




4 


6 




3 


5 




2 


J1054.7-f245 


5 


53 




48 


27 




15 


11 




4 


Jlf03.7-f245 


2 


7 




5 


8 




7 


3 




3 


Jfl22.9-lf36 


f 


8 




3 


10 




4 


13 




4 


Jfl38.2-lf33 


f 


7 




6 


8 




3 


5 







J12f6.8-120f 


2 


32 




9 


23 




9 


11 




9 


J1227.9-ff38 


2 


f4 




13 


28 




11 


5 




5 


Jf354.2-1230 


2 


f5 




12 


12 




4 


11 




9 


All 




f80 


129 


(72%) 


132 


62 


(47%) 


81 


46 


(57%) 



Table 2. Number of galaxies in the full optical field of view con- 
firmed at low or high redshift according to their assumed low or 
high redshift photometric identification. Rows give the candidate 
sub-sample designation, columns give the resulting spectroscopic 
outcome. The percentages in brackets relate to the fraction of the 
candidate sub-sample resulting in either a high redshift or un- 
known classification. The low redshift objects are split into those 
drawn from areas with and without nIR data. The different sub- 
samples are described in section [2.21 



Sample 



Low z objects Unknown 

vith nIR data no nIR data (Low S/N) 



Pri 1 


54 (42%) 


8 


Pri 2 


8 (13%) 


8 


low z 


2 (4%) 


30 


other 


6 (5%) 


45 



10 



17 



57 (44%) 
38 (61%) 
14 (30%) 
45 (40%) 



~400 kms^^, depending upon the location of the features. 
If the continuum break or Lya emission line was compar- 
atively weak, particularly if these features were located in 
the spectral regions containing strong skyline residuals Um- 
iting the accuracy of the redshift determination to ^ 1000 
kms~ , it was classified as grade B. We visually compared 
our sp ectra to the pu b lished one and two-dimensional spec- 
tra of IVanzella et al.l (|2009|) and found that the quality of 
our grade A and B spectra are broadly similar to that of 
similarly classified spectra in their work, even though the 
details of the classification schemes are not identical. 

Over the ten fields, 54 sources were classified as grade 
A 2 > 4.6 galaxies, 36 of which were clear Lya line emitters 
and 18 showed only a Lyman break at the signal-to-noise of 
our data. A further 16 candidates at z > 4.6, but with grade 
B spectra were confirmed, two Lya emitters with weak lines 
and 14 continuum break galaxies. The redshifts of most of 
these placed Lya in spectral regions of comparatively liigh 
sky noise caused by the strong atmospheric lines. Table [5] 
lists the number of high and low redshift galaxies which 
were confirmed within the different candidate sub-samples. 

As can be seen in table (2] the vast majority of galax- 



ies confirmed to be at 2: ~ 5 were drawn from the highest 
priority sample. However, a small number of galaxies which 
had originally been fiagged as likely low redshift interlop- 
ers were found to be at high redshift. Similarly some ob- 
jects fiagged as likely high redshift galaxies were confirmed 
to be at low redshift, highlighting the limitations of purely 
photometric surveys. An accurate redshift for some of the 
confirmed low redshift objects could not be measured due 
to low signal-to-noise or featureless continuum, however a 
2 ~ 5 LBG identification could be ruled out, usually because 
of detectable continuum blue-ward of the expected Lyman 
break at 2 ~ 4.6. 

A further six LBGs were identified in the 'other' sample 
with photometry outside of the selection criteria. Three of 
these have a V-band fiux 0.1-0.2 magnitudes brighter than 
the V-limit (but fulfil the / and R — I selection). Two were 
fainter than the /-band flux limit by 0.2 and 0.3 magnitudes. 
The sixth galaxy had a Rab — Iab = 1.03, bluer than the 
selection limit. This can be explained by the very strong Lya 
emission line detected at 2 = 4.9 from this object. At this 
redshift the line emission falls within the boundaries of the 
_R-band filter, reducing the observed flux diffe rence between 
the R and /-bands (see lStanwav et al.ll2008bl ). 



3.1 Example spectra 

Figure [5] shows an example of a 2 = 5.5 galaxy. At the cen- 
tre of the 2D spectrum there is clear evidence of continuum 
emission with a break across an emission line, a signature of 
Lyman a at high redshift. The continuum break and emis- 
sion line can also be clearly seen in the one-dimensional spec- 
trum. The relatively noisy vertical features in the 2D spec- 
trum are residuals from the subtraction of OH atmospheric 
sky emission. Small-scale apparent continuum features seen 
in the spectrum can all be accounted for by the sky lines, 
atmospheric absorption or noise. 

Figures [3] and |4] sliow examples of low redshift contam- 
inants: an intermediate redshift emission line galaxy and a 
cool star. Several strong emission lines are seen in figure [3] 
as well as detected continuum at the blue end of the spec- 
trum. The line emission was identified as H/3A4861, and the 
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Figure 2. Example of a high rcdshift galaxy at z = 5.5. The 
top panel shows the reduced two dimensional spectrum where 
the horizontal direction is wavelength and the vertical direction 
is the spatial distance along the slit. The bottom panel shows a 
ID extraction from the flux calibrated spectrum and the typical 
atmospheric emission at these wavelengths. 
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Figure 3. Example of a low z emission line galaxy with H/3 and 
the [OIII]AA4959,5007 doublet at z = 0.36. Format as described 
in figure [2] 
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Figure 4. Examp le of an M-type dwarf star (also see 
IStanwav et al.ll2008d V Format as described in figure [2] 










[OIII]AA4959,5007 doublet placing the galaxy at a redshift 
of z = 0.36. 

The stellar spectrum in figure|4]shows a gradual increase 
of continuum into the red, no break, and strong molecular 
absorption bands at 7700A, 8500A and 9300A. At fainter 
magnitudes spectra like this may look like a continuum 
source with a strong break so extra care must be taken when 
identifying Lyman break galaxies without line emission. It 
is more likely that they will appear as a noisy continuum 
source which is comparatively strong at the red end of the 
spectral range, gradually becoming fainter with decreasing 
wavelength. In this case the source would be spectroscopi- 
cally classified as 'unknown' 

Further examples of two-dimensional spectra of line 
emitting and break-only galaxies can be seen in figures [5] 
and IS] 



4 REDSHIFT DISTRIBUTION OF 
CONFIRMED Z ~ 5 GALAXIES 

The redshift distribution of all the spectroscopically- 
confirmed z ~ 5 ERGS sources is shown in figure [T] The 
solid line represents the predicted redshift distribution of the 
candidates given the photometric selection criteria, depth of 
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Figure 5. Two dimensional spectra of a selection of Lya emit- 
ters at increasing redshift with the Lya emission marked by the 
circles. The vertical noise features caused by sky line residuals 
are clearly seen but do not appear to significantly effect the de- 
tection of emission lines. As these spectra have not been rectified 
for individual spectral distortions, the lower Lya redshift scale is 
approximate. 



images and the luminosity function of 2 ~ 5 LBGs as de- 
termined in paper I. The observed redshift distributions are 
shown by the dotted and dashed lines. The grade A galax- 
ies are shown by the dashed line and the whole sample of 
grade A and B objects is shown by the dotted line. The 
general shape of the observed distributions agrees well with 
the prediction, indicating that the parameters of the simula- 
tions (intrinsic colours of the galaxies, shape of the luminos- 
ity function etc.) were reasonable. The atmospheric emis- 
sion, dominated by OH emission lines, at the wavelengths 
of Lya for each redshift is plotted below the redshift distri- 
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Figure 6. Two dimensional spectra of a selection of Lyman break 
galaxies at increasing redshift with the location of the break 
marked by the black lines. Given the faint nature of the objects it 
can be difficult to clearly display the continuum emission in this 
format. The vertical noise features are caused by sky line residu- 
als. As the spectra have not been rectified for individual spectral 
distortions, the lower Lya redshift scale is approximate. 
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butions. Identifying faint liigli redsliift galaxies, especially 
from breaks alone, in regions of liigh noise caused by sky 
line residuals, can be very difficult. For the strong Lya line 
emission and continuum break galaxies (grade A objects) 
this noise does not have a significant effect, but for fainter 
sources areas of high noise can seriously hinder any spec- 
troscopic confirmation. This is reflected in figure [7| where 
the majority of objects with grade B identifications are in 
areas of high sky signal (and hence large residuals after sky 
subtraction) with a lower redshift precision (~1000 kms~^). 
Figure |8] shows the redshift distribution of the sam- 
ple divided into Lya line emitting galaxies and continuum- 
only Lyman break galaxies. Although more line emitters 
were confirmed, the two distributions span the same red- 
shift range. The presence of line emission does not signifi- 
cantly bias the sample completeness as a function of redshift. 
Of the 54 grade A spectra, 36 (67%) had Lya line emis- 
sion. The high fraction of confirmed line emitters compared 
to break-only galaxies is most likely caused by the limited 
signal-to-noise of the spectroscopy. The redshift of compar- 
atively strong line emitters is much easier to determine than 
the location of a continuum break. When the whole sample 
of grade A and B spectra of 70 galaxies is considered 38 
(54%) confirmed high redshift galaxies have Lya line emis- 
sion, resulting in a more balanced detection rate for both 
populations. We stress that a grade B rating does not indi- 
cate uncertainty in the high redshift identification, only in 
the precision of the derived redshift. 



Figure 7. Redshift distribution of spectroscopically confirmed 
galaxies. The dotted line represents all galaxies with rcdshifts 
z > 4.5, the dashed line shows only grade A galaxies with secure 
redshifts and the solid line illustrates the shape of the predicted 
redshift distribution (scaled to match the peak number count of 
the observed distribution). Also included, below the distributions, 
is the atmospheric emission spectrum at the redshift of Lyo. 
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5 REST-FRAME UV CONTINUUM SLOPE 

In the local universe, the spectral slope in the UV, /3 
(where Fx (x A'^), is sensitive to the amount of redden- 
ing, the fraction of the non-ionising continuum that is re- 
process ed into infrared emiss ion, and the metallicity of the 
galaxy (|Heckman et al.lll998l ). There are good correlations, 
with similar significance, between the spectral slope and 
log L/_r/L[7v, the strength of the relatively low ionisation 
metal lines from the ISM, and the average metallicity. In 



Figure 8. Redshift distributions of Lyo line emitting LBGs (dark 
grey) and continuum only LBGs (plotted above the line emitters 
in light grey). Both grade A and B LBGs are included. Also in- 
cluded, below the distributions, is the sky emission spectrum at 
the redshift of Lyo. 
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addition /3 is relatively insensitive to the age of the galaxy 
within the first 100 Myrs, lying between —2.6 and —2.0 
(JLeitherer fc Heckmanlll995l ^. These correlations reflect the 
complex relationships between the population of stars that 
emit significantly in the non-ionising UV (dominated typ- 
ically by B-stars but is a strong function of age for the 
youngest populations, e.g.. iMeurer et al.l 119971 ). the over- 
all extinction and reddening, and the metallicity of the gas 
(which infiuences the gas to dust ratio, the gas temperature 
and ionisation state, and the metal line absorption strength). 
These relationships can be used to investigate the likely 
properties of our distant galaxy sample. We can also test 
the appropriateness of these local relationships by making a 
comparison with z~3 LBGs. 



5.1 UV slope from spectroscopy 

Two simple average spectra were created, one for the Lya 
emitting sources and one for the break-only LBGs, by com- 
bining the sources with the best signal-to-noise continuum 
detections just long- ward of 1216A in the rest-frame us- 
ing redshifts determined by the location of the continuum 
break and line emission. A one-dimensional spectrum was 
extracted from each two-dimensional flux calibrated spec- 
trum, converted to a rest-frame wavelength scale and com- 
bined with similar spectra of other galaxies with or with- 
out line emission. Spectra of ten sources were combined to 
form the average line-emitter spectrum and thirteen for the 
average break-only spectrum. The difference in the num- 
ber of spectra used for each average affects the signal-to- 
noise, but not the normalisation or shape of the final spec- 
tra. The two sets of sources had comparable redshift distri- 
butions, average /—band magnitudes and signal-to-noise at 
~ 1250A in their rest-frame spectra. All spectra were drawn 
from the same sets of masks and identically reduced with the 
same calibrations; there was no way of determining a priori 
whether an object would be a line-emitter or a break. While 
the line-emitters were identified from the presence of both a 
line and break, the identification of the break-only sources 
required a clear continuum break across 1216A, rather than 
a gradual fading towards the blue and so could potentially 
bias the highest signal-to-noise break-only sample used here 
to those galaxies with the bluest breaks. 

A comparison of the rest-frame UV slope for the two 
LBG populations is shown in flgure[51 Although this compar- 
ison indicates that break-only galaxies have a slightly redder 
continuum than the line emitters (and show a difference in 
the opposite sense to the potential bias noted above), this 
could be due to the Ly« absorption or broad dust features 
rather than the intrinsic slope of the continuum. Without 
near-infrared spectroscopy, sampling longer wavelengths, re- 
liable measurements of the UV slope cannot be made with 
this data. However, at the resolution and wavelength cov- 
erage of the 2: ~ 5 spectra presented here, both stacks are 
consistent with sim ilar stacks of l ine-em itting and break- 
only z - 3 LBGs in lShaplev et all (|2003l ). 



5.2 UV slope from photometry 

Prev ious studies have used I—z colours fe.q. lPentericci et al.l 
120071 ). in the absence of suitable spectroscopy, to measure 
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Figure 9. Average spectra of line emitting (solid line) and break- 
only galaxies (dashed line) smoothed using a boxcar algorithm 
with a width of 5A The averages were made from spectra of the 
10 grade A line emitters and 13 grade A break galaxies with the 
highest S/N continuum detections just longward of rest-frame 
1216 A, the two sets of galaxies with similar redshift distributions 
and average /—band magnitudes. 



the slope of the rest-frame UV continuum of LBG galaxies 
at z = 3 — 5. At 2 -^ 5, the ERGS z-band imaging (and spec- 
troscopy) does not provide a wavelength baseline (rest-frame 
1200A to 1450A) sufficient for estimating the continuum 
slope reliably. Instead, the longer wavelength Ks-h&nA data 
is used in combination with the /—band, resulting in a rest- 
frame wavelength range of ~ 1300A to --^ 3500A. Although 
the near-infrared data in this study is not deep enough to 
detect individual high redshift galaxies, using the technique 
of stacking many objects together enables the data to reach 
fainter magnitudes. 

Sections of the /sTg-band images, centred on the objects 
of interest, were stacked together using imcombine within 
the IRAF package. The mean stack was created with per- 
centile clipping and was then used to measure the average 
Ks-ha,nd flux. Two examples of the LBG stacks are shown 
in figure [TOl the Lya line emitting and break-only samples. 
At 2 ~ 5 the Lyman break falls into the /-band which af- 
fects the measurement of the continuum level. To correct 
for this effect and obtain a colour which represents only the 
continuuiu, the /-band luagnitude of each galaxy contribut- 
ing to the stacks was corrected for the redshift dependent 
continuum break effect and Lya line emission. This gives an 
average Iab — Kab colour for the whole sample which can be 
transformed into a measure of the UV to blue optical spec- 
tral slope. Table [3] summarizes the stacking of the Ks-hanA 
imaging data of spectroscopically confirmed 2 ~ 5 galax- 
ies across the ten survey fields. We found that the break- 
only galaxies are significantly redder than the line emitters, 
Iab — Kab = 0.8 compared with Iab — Kab < 0.0, the 2a 
colour limit. 

Given the degeneracy between stellar age, reddening 
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Table 3. Stacked ly^g — K^g colour of confirmed high redshift 
galaxies, break-only galaxies and Lya line emitting galaxies. The 
individual /-band magnitudes of the galaxies have been corrected 
for line emission and the redshift dependent location of the Ly- 
man break within the /-band filter. The rest-frame UV slope is 
described by /3 assuming F^ <x A''. Where a limit is shown, the 
2cr depth of the stacked Ks-hand was used. 



Sample 


No. gals 


j-mean 


Iab — Kab 


P 


All breaks 
All lines 


20 
22 


25.56 
25.62 


0.82 ±0.23 
<0.02 


-1.3 
<-2.0 


All z ^ b galaxies 


42 


25.59 


0.27 ±0.25 


-1.8 










Figure 10. Imaging stacks of the /C-band images of break-only 
galaxies (left), Lya emitters (right). Both images have the same 
smoothing and scaling applied. 



and metallicity it is difficult to determine the primary source 
of the observed colour difference between the two subsam- 
ples. Nevertheless, we compare the relative effects of stellar 
age and redde ning to the Iab — Ka b colour. 

Using the lCalzetti et al.l (|2000l ) dust models, the colour 
difference between the two subsamples of at least Iab — 
Kab ~ 0.8 can be produced by a small difference in the av- 
erage redenning between them, with the break-only galaxies 
being reddened by an extra E(B-V)~ 0.2 relative to the Lya 
emitters. The average colours of both samples indicate a low 
level of dust content, with a typical E(B-V) of more than a 
few tenths, otherwise the intrinsic colours would be too blue 
for even very young models with a standard IMF. The levels 
of reddening are consistent with those found in photometric 
studies of z ~ 5 LBGs (e.g. IVerma et al.ll2007l) and in spec- 
troscopic studies of z ~ 3 LBGs {e.g. IShaplev et al.ll2003l . 
finds E(B-V)=0.15 and 0.1 for break-only and line emitting 
galaxies at 2 ~ 3). 

The colour difference could also indicate an average stel- 
lar population age difference between the two subsamples, 
though any age difference is constrained by the age of the 
Universe a.t z ^ 5 (about 1.1 Gyr) and, more stringently, 
with the range of ages found in previous photometric stud- 
ies of individual z ~ 5 galaxies. For example, the study by 
IVerma et al.l (|2007l ) used detailed photometry out to 8/im to 
determine that the majority (~ 70%) of such sources have 
young (< 100 Myr) dominant stellar populations, with only 
a minority older than this, typically up to 2-300 Myr. 

Using composite stellar population models from 



iMarastonI (|2005l ) that bracket the plausible range of star 
formation and metallicity for the galaxies, we explored how 
age differences could affect the observed Iab — Kab colours. 
For an exponentially declining star formation rate, with an 
e-folding time of 30 Myr and a metallicity of 0.05Zq, a 
colour difference of Iab — Kab ~ 1 could be obtained if 
the youngest subsample was no older than a few 10s of Myr 
and the oldest ^ lOOMyr older. If the youngest population 
was itself on average > 100 Myr old, the age difference in- 
creases to 200 Myr or more. For a model where the e-folding 
time was 100 Myr and the metallicity 0.5 Zq, the older sub- 
sample had to be typically 200 Myr older than the younger 
for almost any plausible age of the younger population out 
to 500Myr. For an unreddened Iab — Kab ~ 0.8, a source 
would need to be 200 or 400Myr old for an e-folding time 
of 30 and lOOMyr respectively. An unreddened colour limit 
of Iab - Kab < results in ages of < 100 and < 200 Myr 
for the same models. Associating the line emitters with the 
youngest populations and the break-only galaxies with the 
older, the ages of the old er population s are d ifficult to recon- 
cile with the results from lVerma et al.l(|2007h . given that half 
of the spectroscopically-confirmed sample are break-only 
systems. Clearly it is possible that a combination of differ- 
ences in both stellar population age and reddening could ex- 
plain the observed colour difference between th e subsamples, 
but as suming the age constra ints indica t ed bylVerma et al.l 
(J2003) and the colours of the IMarastonI (|2005h models, the 
colour difference is more simply achieved through a variation 
in redenning. 

A difference in the intrinsic rest-frame UV continuum 
slope for LBGs has been seen in other photometric and 
spectroscopic studies at ^ = 3 — 5. At z = 3, UV slopes 
oi P — —0.7 and /? — —1.1, for break and line emit- 
tin g galaxies resp e ctively , are measur ed from spectra l stacks 
in Ishaplev et al] l|200a) . Similarly IVanzeUa etU] (|2009l ) 
presents a z ~ 4 spectral stack for the two galaxy pop- 
ulations which is c onsistent with the UV slopes measured 
photometrically in iPentericci et al.l (|2007l ) (/3 = —1.7 for 
break galaxies and /3 = —2.0 for line emitters). The results 
presented here follow similar trends with a photometrically 
measured slope of -1.3 for break-only galaxies and <-2 for 
Lya emitting galaxies. 

Given that no rest-frame UV colour criteria were used in 
our observational selection {e.g. no Iab — zab colour cut), 
the observed targets were not biased to objects with blue 
continuum. However, the majority of confirmed LBGs were 
observed to have relatively blue, fiat (/3 ~ —1.8) rest-frame 
UV continuum. Only two galaxies were detected individ- 
ually in the near-IR imaging suggesting a redder contin- 
uum. As the objects observed are bright in the UV, it is 
to be expected that they have blue continuum dominated 
by young stellar populations. If a sample containing objects 
with fainter UV luminosities were probed, it would likely 
contain a higher proportion of older and perhaps more red- 
dened galaxies. 

5.3 The Metallicity of z~5 LBGs 

iHeckman et al.l ([1993) found correlations between the spec- 
tral slope and log ^ir/I^uv, the strength of the relatively 
low ionisation metal lines from the ISM, and the aver- 
age metallicity for nearby galaxies. In Figure 1111 we plot 
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consistent with forming at z ~ 5. The growth in the number 
density and in the typical mass of a LBG over this period are 
also roughly a factor of 10. This may suggest a direct evolu- 
tionary connection between the two populations (cf. paper I, 
IVerma et alll2007l : IShaplev et allboosf : ISteidel et aLlll999l) 
Moreo ver, in an analysis of z ~ 3 LBGs, iMannucci et al.1 
(120091) find that galax ies with masses si milar to those ex- 
pected for this sample (|Verma et al.ll2007l ) have metallicities 
consistent with our estimates using the measured z ~ 3 /3 
and the local relationship between /3 and metallicity. For 
such systems, IMannucci et al.1 ((20091) have estimated gas 
fractions of 80% or more. For the oldest s ystems at z -^ 3 , 
which are also the most massive (see also lErb et al.ll2006l ). 
the gas fractions are estimated to be much less (~20%). 
Therefore, although beyond the scope of the present paper, 
it is possible to relate the more massive, older, more metal- 
rich populations of LBGs at z ~ 3 to the population we 
observe at z ~ 5. 



Figure 11. Plot of the UV spectral slope, /3, versus the metal- 
licity. The black d ots show the spectral s lope and metallicity of 
local galaxies from iHeckman et al.l 1119981 ). the dot-dashed verti- 
cal lines indicate the spectral slope of the galaxies with breaks 
only in our sample (red) and for the whole sample with spec- 
troscopic redshifts (blue), and the dashed vertical lines indicate 
the spectral slope of the galaxies with Lyo mostly in absorption 
(red) and galaxies with L ya predominately in emission (blue) 
from lShaplev et al.l 1 120031 ). The red horizontal region shows the 
area encompassing the metallicity estimates from the literature 
(see text for details). 



the spec tral slope, ff, versus t he metallicity, 12+log[0/H], 
from the lHeckman et al.l (| 19981 ) sample of local galaxies and 
their best fit to these data. Overplotted are the estimates 
of the spectral slope of z ~ 3 LBGs from IShaplev et al.l 
((2003) and t he range of m etal abundances determined for 
z~2, LBGs JPettini et"ai][2b01: : Ncsvadba ct al. 2006 , 200^ ; 
Mannucci et al.ll2009l ). Two values of the spectral slope from 
Shaplev et al.l (|2003l ) were used, a steeper slope for the Lya 
emission line galaxies and a redder slope for the galaxies 
where Lya is in absorption. 

We see that, in general, /3 for the 2: ~ 3 LBGs pre- 
dicts a higher metallicity than has been observed. The num- 
ber of good metallicity estimates is low and the comparison 
may suffer from biases (such as slope estimates weighted 
towards the brightest galax;ies, which have populations that 
are more evolved, or the emission line estimates are weighted 
towards galaxies of lower metallicities and dust content on 
average due to the influence of metallicity on the gas tem- 
peratu re, effecting the line s trengths of rest frame optical 
lines). iHeckman et al.l (| 19981 ) found that galaxies brighter 
in the optical and with higher (approximate) bolometric lu- 
minosities have bluer spectral slopes. Realising that in this 
work we are discussing the brightest LBGs at z~5, the spec- 
tral slopes suggest that the z ~ 5 LBGs are roughly a factor 
of three lower in metallicity than the 2-^3 LBGs and are 
approximately 1/5-1/lOth of solar metallicity but this esti- 
mate may be biased towards galaxies of higher metallicities. 

About 10-20% of the LBGs at z ~ 3 have stellar ages 



6 LYMAN a LINE PROPERTIES 

The Lyman a line is the only emission or absorption line 
straightforwardly detectable in the observed optical in a rea- 
sonable time frame (four hours using a 8m telescope) in 
the spectra of LBGs at z ~ 5. A small number of observa- 
tions of the continuum of bright or lensed z ~ 5 galaxies 
have reached sufficient depths to probe interste llar absorp- 
tion features in the rest-frame UV co ntinuum (|Ando et al.l 
|2004 I2OO7I : iDow-Hvgelund et al.|[2005l ) . 



6.1 Lyman a at z ~ 5 

The Lya line flux of each galaxy was measured within 
a square aperture centred on the line emission in the 
two dimensional spectrum. The size of the aperture was 
tuned to the apparent spatial and spectral FWHM of each 
galaxy. The flux was integrated over the wavelength range 
within the aperture, yielding a final measurement in units 
of ergscm"^ s"'^. Line fluxes ranged from 2.6 x lO"'^^ to 
7.0 X lO^^^ergscm^^ s^^ (luminosities between lO^'^ and 
10*'' erg s^^) with a typical 3o- line flux limit of 2.9 x 
10^^* ergscm^^ s^^ for an unresolved emission line. At the 
low redshift end of the sample there is a lack of strong emis- 
sion lines as a result of the line flux contributing to the R- 
band magnitude, reducing the R — I colour and causing the 
object to drop out of the original colour selection. From the 
z > 5 Lya equivalent width distribution potentially 10-30% 
of line emitters could be missed at z < 5 give n the influence 
of Lyg emission on the Rab — Iab colour (jStanwav et al.l 
l2008bl ) 



6.2 Equivalent Width of Lya 

As the observed galaxies are faint (26.3 > Iab > 25.0), 
the continuum emission detected spectroscopically has low 
signal-to-noise. As such, equivalent width measurements (ef- 
fectively the ratio of the Lya line flux to continuum flux den- 
sity) given here use the broadband photometry of individual 
galaxies to estimate their continuum flux. Measured line flux 
was removed from the 2" aperture 7-band magnitude to de- 
termine each galaxy's continuum flux level, accounting for a 
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Figure 12. Distribution of tiio rest-frame Lyo equivalent widtii 
for line emitters (solid line) and when 3 a upper limits from 
spectroscopically-confirmed continuum break galaxies are in- 
cluded as described in the text (dashed line). Continuum flux 
levels were measured from broadband photometry. 



Figure 13. Estimated distribution of rest-frame Lya equivalent 
width upper limits for observed sources without distinguishable 
features in their spectra. This distribution was determined by 
convolving 3 cr upper limits measured from the spectra as a func- 
tion of rodshift with the redshift distribution calculated for the 
sample, as described in the text. 



redshift dependent correction due to IGM absorption short- 
ward of Lya. A limit on the Lya equivalent width of the 
break-only galaxies was measured by introducing an artifi- 
cial hne feature with a flux equivalent to the 3 a variations 
in the local background at the break wavelength. This sim- 
ulated emission line was placed at the continuum break of 
the spectrum and used to estimate an upper limit on the ob- 
served equivalent width limit which accounted for the noise 
properties of each spectrum. The derived rest-frame equiv- 
alent widths are shown in figure 1121 The strongest line de- 
tected has a rest-frame equivalent width Wo = 112 ± 18 A, 
the weakest detected line had Wo = 6 ± 2 A. 

While figure [12] presents the distribution of Lya line 
strengths for sources with spectroscopically confirmed red- 
shifts, this is an incomplete picture. To develop a more com- 
plete understanding of the equivalent width distribution at 
z ~ 5, it is necessary to consider those sources which were 
spectroscopically observed, but did not yield a reliable red- 
shift determination or recognisable spectral features. Unsur- 
prisingly, these sources outnumber those with spectroscopic 
confirmations and must comprise a mixture of faint galax- 
ies, sources at intermediate brightness with weak or absent 
emission lines, and interlopers contaminating the photomet- 
ric sample. 

In Paper I we determined that ~80% of our refined 
photometric sample is likely to lie at z > 5. Considering a 
similar subset of sources with reliable multiband photome- 
try, selected as strong high redshift candidates and placed 
on the spectroscopic masks, we identify 57 sources for which 
no redshift was determined, of which ^ 45 are expected to 
be genuine high redshift galaxies, but whose emission lines 
may lie anywhere in the observed spectrum. In the follow- 
ing we estimate the contribution of the upper limits to any 



Lya emission from these sources to the complete equivalent 
width distribution. 

We calculate a probability distribution of upper lim- 
its for this sample as follows. For each source in the sam- 
ple, we determine a 3 cr upper limit on the flux of an unde- 
tected emission line as a function of redshift, accounting for 
wavelength-dependent noise in the spectrum. This is con- 
verted to a limit on the rest-frame Lya equivalent width as 
a function of redshift, using the 7-band magnitude of the 
source to determine a continuum level (and accounting for 
IGM absorption shortwards of the proposed emission line, 
again as a function of redshift). We weight the probabil- 
ity of the line occurring at each redshift by the calculated 
redshift distribution function of this sample (figure[7]), yield- 
ing a probability distribution for the 3 a upper limit on the 
equivalent width for each source. The distribution is then 
scaled by 80% to produce an estimated contribution to the 
total sample equivalent width distribution. The resultant 
distribution of line upper limits is shown in figure 1131 The 
typical 3(7 upper limit on a line detection is 5 ± 2 A, with 
the small amplitude of the tail to higher limits suggesting 
that our sample is highly complete for line emitters down to 
our flux limit. 

Any interpretation of the Lya line is challenging given 
its complex and ambiguous nature. The observed line profile 
is affected by resonant scattering within the galaxy by in- 
terstellar HI. This introduces influences from the geometry, 
kinematics and dust content of large-scale outflows. 

How does the z ~ 5 equivalent width dis tribution com- 
pare t o that at lower redshifts? At 2; ~ 3, IShaplev et al.l 
([2003) determined the equivalent width distribution of Lya 
for nearly 1000 galaxies. They flnd that 25% of 2; ~ 3 Lyman 
break galaxies have a Lya equivalent width greater than 
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20A. Of galaxies in our spectroscopically-confirined sample 
at z ~ 5, 55% have Wo > 20A. However, when we con- 
sider the larger sample of spectroscopically-observed high 
redshift candidate sources and the estimated contribution 
of unconfirmed galaxies in figure [13] is taken into account, 
this fraction falls to 22 ± 4%. 

Previous work at lower redshift including lShaplev et alJ 

(|2003l ) similarly did not include a correction for unconfirmed 
sources. While the confirmed fraction at z ~ 3 (63%) in that 
work is significantly larger than that in this sample, and 
included sources with Lya seen in absorption, it is likely that 
an essentially unknown fraction of the unconfirmed sources 
would contribute to the Lya equivalent width distribution 
in the same manner as those discussed above. Consequently, 
the quoted fraction of sources with high equivalent widths 
at z ~ 3 is an upper li mit, but nevertheless , even if all 
unconfirmed sources in the lShaplev et al.l (|2003l ) sample are 
at « ~ 3, the high equivalent width fraction drops to ^ 18%, 
consistent with our z ^ 5 statistics. 

Interestingly, this survey is not alone in finding a large 
fraction of confirmed high redshift galaxies with a large Lya 
equivalent width. A long tail of source s with high equivalent 
width s is seen bv lDawson et al. (12007 ) at z = 4.5. IHu et al.l 
l|2004h aX z = 5.7 and IStanwav et all l|2007h at z ~ 6, al- 
though these studies disc uss only sources with spectr oscopic 
confirmation. In contrast lDow-Hvgelund et al.l (|2007l ) do not 
observe such a tail based on 6 line emitting galaxies at z ~ 6. 
At present, all of these studies suffer from low spectroscopic 
completeness and small number statistics. The current study 
is the first at high (2 > 5) redshift with sufficient number 
counts and understanding of selection functions to make a 
meaningful comparison of the equivalent width distribution 
to that seen at lower redshifts. 

The largest equivalent widths seen at z ~ 3 are almost 
double those at 2: ~ 5. Given their rare occurrence, their 
absence in the higher redshi ft sample may simply be an effect 
of small number statistics. IShaplev et al.l (|2003. ') also found 
that the UV continuum slope became bluer as the equivalent 
width of Lya increased in stacked samples of LBGs at z ~ 3. 
Figure [TJ] shows the distribution of Iab — zab colour, a 
rough indicator of the UV continuum slope once corrected 
for line emission, with the rest-frame equivalent width of 
Lya. Although a comparably clear trend is not seen, sources 
with equivalent widths of more than 30A appear to have 
bluer colours than those with lower equivalent widths. There 
are more z-band detections resulting in redder colours for 
sources with Wo <30A. However, given the small number of 
sources this is not highly statistically significant. 

The resolution of the spectroscopy was such that only 
the strongest Lya lines were barely resolved, so little can be 
said about the distribution of line widths and line shapes in 
the sample; higher resolution spectroscopy is required. 



7 MORPHOLOGY OF Z ~ 5 LBGS 

HST/ACS imaging using the F814W filter (/-band, 
iDesai et a l. 2007) were used to investiga te the morphology of 
the confirmed high redshift LBGs fsee lDouglas et al.ll2009l . 
for details). Of the spectroscopically confirmed LBGs, 20 in- 
dividual sources have deep HST imaging. Within this sam- 
ple, 14 LBGs have multiple components, four show irreg- 
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Figure 14. Rcst-framo equivalent width of Lya as a function of 
Iab — ZAB colour. Diamonds represent equivalent width limits 
of the break-only galaxies as described in the text. The colour 
has been corrected for line emission and IGM absorption effects. 
Arrows indicate a non-detection in the 2-band and hence a colour 
limit. 



ular morphology, one is compact and unresolved and the 
final object is compact but surrounded by diffuse emission. 
Examples of LBGs with multiple components or irregular 
morphology are shown in figure [151 The morphology of indi- 
vidual spectroscopically confirmed high redshift galaxies is 
similar to that of the larger photome trically selected sam - 
ple of z ~ 5 Lyman break galaxies (JDouglas et al.ll2009l \ 
A similar mean half-light radius of 0.15" (or 0.9kpc) and 
distribution of sizes was observed for the two samples, and 
are in agreernent with other high resolu tion studies {e.g. 
[Bremer et al] [20041 : [Bouwens etahl [2006h . The number of 
components identified in the multiple systems ranged from 
two to four, however the systems may contain additional 
fainter components that are not detectable above the sur- 
face brightness limits of the images. The mean distance be- 
tween components is 0.84" (5.3 kpc) with distances observed 
between 0.2" and 2.5" (1 kpc and 16 kpc) compared to a see- 
ing disk of typically 0.7" in the groundbased imaging. 

Two systems are resolved in the ground-based imaging 
with component separations of 0.8" and 1.0", however fur- 
ther components have been identified in the deeper, high 
resolution HST images. Within each system, the previously 
identified components have similar optical colours with no 
detections in the nIR imaging. Both systems are marginally 
resolved in the spectroscopy showing two sources. One sys- 
tem has two faint Lya emission lines and the other has a line 
emitter and a break-only galaxy. The component redshifts 
measured from these features are within 1000 km s^ ^ , which, 
although high, does not rule out a connected system as red- 
shifts measured using Lya emission or the Lyman break are 
not necessarily equal to the systematic redshift of the galaxy. 
Velocities of up to 1000 km s~^ between Lya emission and 
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Figure 15. Eight examples of spectroseopically eonfirmed galax- 
ies with multiple components. Each box is 5" X 5" and the im- 
ages were taken using the HST/ACS camera and the F814W fil- 
ter. Redshifts are indicated in the bottom right-hand corners 



interstellar absorption lines have been observed in z ~ 3 
LBGs. 

A higher incidence of complex morphology in the 
spectroscopically-confirmed sample is observed compared to 
the larger photometric sample where ~30% of galaxy candi- 
dates had more than one component or irregular morphol- 
ogy. This effect has not b een observed in other spect roscop- 
ically confirmed samples (jVanzella et al.ll2005l . 12009 ). How- 
ever, given the small number of spectroseopically confirmed 
zr;5 LBGs and differences in selection and image depths, 
these numbers are statistically consistent with other s tudies 
IVanzella et al.ll2005l . [2009l : IConselice fc Arnoldll2009l '). 



7.1 The nature of the multiple components 

The relatively high fraction of distant LBGs that have mul- 
tiple components may tell us important information about 
their nature of their star-formation and its triggering mech- 
anism. The high fraction of pairs or triplets suggests that 
the mergers and instabilities that form massive regions of 
intense star-formation within larger underlying structures 
may play an important role in z ~ 5 LBGs. The observed 
faint structures surrounding these sources might be indica- 
tive of mergers, but could also be structure within the light 
distribution of a single larger object. There is currently no 
direct evidence to favour either scenario. 

An appropriate way to understand the nature of Z'-^S 
LBGs may be through analogy with a local population. 
lOverzier et al.l (|2008l ) have studied nearby analogues of 
LBGs - Lyman Break galaxy analogues or LBAs. These 
LBG analogues were chosen to have properties similar to 
z ~ 3 LBGs, namely, Lfuv > W^°'^ Lq and Ifuv > 10^ 
Lq kpc~^. Thus they mimic the high UV luminosities and 
intensities observed in z ~ 3 LBGs. Overall, the LBAs ap- 
pear to be similar to the z ~ 5 LBGs, perhaps more so 
than the z ~ 3 LBGs they were choosen to mimic. Their 
total stellar masses are between 10^ and 10^" ^Qy close to 
the typical masses of bright (brighter than L*) z ~ 5 LBGs 
IjVerma et al.ll2007l ). and low metallicities and similar spec- 
tral slopes ( §[531 and §E21 ^\EM- Their UV luminosities 
are typically a factor of a few to 10 times lower (APa for 
A=1900A) than our sample. Both the small scale clustering 
properties and the internal kinematics of the LBAs suggest 
a merger or interaction ori gin for the intense star- formation 
observed in these systems (jBasu-Zvch et al.ll2009al lbl) . Given 



the apparent similarities, a detailed comparison seems war- 
ranted^ 

lOverzier et al.l l|200a ) find that the clumps of star- 
formation, or SSBs, in the LBAs, have mass surface den- 
sities of lO'^"^ Mq pc~^ over radii of ~100 pc, with typi- 
cally one to more than 10 large SSBs per galaxy. Such sys- 
tems would be unresolved at z~5 in any of our datasets. 
The half-light radii of the galaxies themselves in the UV 
or U-band are typically around one kpc. A comparison of 
this mass surface dens i ty wi th the average surface densi- 
ties of the lVerma et al.l (|2007l ) sample (comparable to LBGs 
studied here) gives a similar average mass surface density 
(~6xl0^ Mq pc~^ over radii of ~1 kpc or the half light 
radu of z ~ 5 LBGs). In the LBAs, about 30-80% of the 
total inferred star-formation occurs in such super-starburst 
(SSB) regions. Thus, the z ^ 5 LBGs could be thought of 
as a "galaxy wide" super-starburst region, with a high den- 
sity of star-formation that must contribute nearly 100% of 
the inferred star-formation rates (which would also explain 
the higher UV luminosities for the z ~ 5 LBGs compared 
to their local analogues). In addition. lOverzier et al.l (|20Q8l ) 
infer ages of the SSBs of about 10-100 Myrs, but formally 
regard these as upper limits due to complex extinction in the 
UV and due to the (unkn o wn) c ontribution of older stars. 
Interestingly, IVerma et al.l ((20071) derived similar ages from 
model fits of the bandband photometry of z ~ 5 LBGs. This 
is all consistent with the picture that the z ~ 5 LBGs are 
possibly a collection of galaxy wide SSBs with a contribution 
to the total star- formation of nearly 100%. 

If the LBGs studied here are analogous to the LBAs, 
with spatial distributions over one to 5 kpc, with typical 
separation between components of around 1 kpc, we would 
only expect to distinguish so me of the separate SS B clumps 
in HST images {e.g., Fig. 7 in lOverzier et al.ll2008h . For the 
spectroscopically-confirmed galaxies with deep HST imag- 
ing, we find 12 that are clearly multiple with an aver- 
age projected separation of 5.3 kpc and range from sepa- 
rations of one to 16 kpc. For the total photometric sam- 
ple with deep HST imaging, 14 out of 95 are multiple sys- 
tems and 16 out of 95 show disturbed mor phology (mod- 
ulo t he contamination which should be low iDouglas et al.l 
[2OO9I). Taken at face value, the high fraction of multiple 
or disturbed systems suggests a bias towards spectroseopi- 
cally confirming such objects. This suggests that the galaxies 
with multiple components of star formation allow for chan- 
nels to be cleared of neutral hydrogen, which then provide 
a long mean free path for the Lya photons to escape. Since 
these galaxies likely drive outflows (jVerma et al.ll2007l ). the 
nature of the ISM may be such to allow the efficient es- 
cap e of Ly g photo ns in analogy with local starburst galax- 
ies (iMarlow c. Mcur er. fc Heckman|[l999l : iMeurer et al.ill993 : 
lAtek et al.li2008. : .Haves et al.ll2007l ). 



8 Z ~ 5 SURFACE DENSITY 

The average observed surface density of the most reliable 
photometric sample of galaxies was 0.93 per arcmin'^ down 
to Iab = 26.3 (see table 2 in paper I), comparab le to the sur- 
face density observed bv lBouwens et al.l ([20071) to the same 
flux limit from space-based imaging surveys. However, there 
is a large variation in the surface density of LBG photo- 
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Table 4. All confirmed high redshift Lyman break galaxies per 
field according to the quality of the spectra and the identifying 
feature. 
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8.1 Spectroscopically confirmed overdense fields 

Out of the ten survey fields, two showed strong spectroscopic 
evidence of large scale structure. In fields J1040. 7-1155 and 
J 1054. 7- 1254 there are clear spikes in the redshift distribu- 
tions at z = 5.15 and z — 5.0 respectively, with six or more 
galaxies lying within a 0.1 redshift bin, which is equivalent 
to a co-moving distance of around 50Mpc. The redshift dis- 
tributions for these two fields are shown figures [T6l and [T9l for 
J 1040. 7- 1155 and J 1054-7. 1245 respectively The red dashed 
line indicates the distribution of grade A galaxies, the blue 
dotted line shows the distribution of all grade A and B galax- 
ies, and the black line shows the night sky emission spectrum 
for the redshift of Lya. The total redshift distribution of all 
confirmed galaxies predicts, at most, 1.3 galaxies per 0.1 
redshift bin for an individual field. 



metric candidates and spectroscopically confirmed objects 
between individual fields. Table U show this variation in the 
number of spectroscopically confirmed objects. As expected, 
the number of confirmed galaxies in each field broadly cor- 
relates with the number of photometric candidates. A factor 
of 8 is observed between the richest and poorest LBG candi- 
date fields and those with the most and least spectroscopic 
identificatio ns, beyond that expecte d from cosmic variance 
predictions JTrenti fc Stiavelli|[2008l '). 

One possible explanation of the observed variation is 
the changing effective area of each field. However, with vari- 
ations in area of less than 10% and no correlation with the 
density of candidates and confirmed LBGs suggests that this 
cannot account for the differences. 

Alternatively, the presence of foreground clusters in 
each field (the original subject of the maj ority of the 
imag ing by the ESO large program EDisCS, ([White et al.l 
[2OO5I)) could be di storting the sourc e plane and affecting 
the probed volume (|Clowe et al.ll200(i ). Given the steep na- 
ture of the faint-end of the z ~ 5 LBG lu minosity function 
iDouglas et al.ll2009l : lBouwens et al.ll2007l ). the number den- 
sity of candidates could increase significantly with the fiux 
boosting effect of lensing. As the area of each field affected 
by strong lensing is small, this is unlikely to be the cause 
of the observed number variation. There are also no objects 
observed with highly distorted or elongated morphologies, 
suggesting strong lensing in the deep images. Although the 
effects of weak lensing cover a larger proportion of the field, 
the subsequent distortion effects are small and do not signif- 
icantly alter the probed volume. There is also no correlation 
seen between the number of candidates and confirmed LBGs 
and the velocity dispersion, and hence mass, of the identified 
foreground clusters. Therefore it is unlikely that the lensing 
effects of foreground structures are causing the observed dif- 
ferences in LBG surface densities. It is far more likely that 
these fields are probing truly overdense areas of the universe. 
The most compelling evidence for this is the variation from 
a smooth redshift distribution seen in the spectroscopic re- 
sults of two of the individual fields and is discussed in the 
next section. 



8.1.1 JIO4O.7-II55 

This field contains the third largest photometric sample of 
high redshift galaxy candidates out of the ten survey fields, 
with 25 candidates in the photometric sample (paper I) and 
11 high redshift spectroscopic confirmations (grade A and 
B). Given the size of the photometric sample, a similarly 
high number of spectroscopic confirmations would be ex- 
pected. In the redshift distribution there is a clear peak of 
objects between 5.11 < z < 5.21, a 6a overdensity when 
compared to the spectroscopic statistics of the survey as 
a whole (with the boundaries of the redshift bins tuned 
to give the largest peak). Figure [161 shows this peak. Only 
line emitters with redshift determinations accurate to within 
300-400 km s~^ populate the redshift spike strengthening its 
identification as a real feature. It cannot be explained as 
an artifact created by the uncertainty in a Lyman break 
overlying skyline residual noise. Figure [T7] shows the spa- 
tial distribution of the galaxies in the redshift spike across 
the field. Although the galaxies are spread across the field 
and are not concentrated in one region, it is still probable 
that they are part of the same structure. The closest pair of 
grade A galaxies in the field are 7.3 Mpc apart, with a fur- 
thest extent of 43.7 Mpc for the spectroscopically confirmed 
grade A members of the cluster, where the distances are co- 
moving in three dimensions. The co-moving dimensions of 
the surveys fields are typically 27 by 27 Mpc in the plane of 
the sky and probing a redshift interval of one, equivalent to 
525 Mpc. 

The spectra of the galaxies constituting this redshift 
spike are typical of the line emitting sample as a whole. Sim- 
ilarly, the morphologies and colours of the galaxies within 
the redshift spike are typical of the sample. 



8.1.2 JIO54.7-I245 

This field contains the most photometric high redshift can- 
didates (45 candidates in paper I), and was expected to 
produce the largest number of spectroscopic confirmations. 
Within the spectroscopic identification of 24 high redshift 
galaxies (grade A and B), there is significant spike at 2 = 5.0 
within Sz = 0.1, containing six times more galaxies than ex- 
pected over such a redshift range (grade A only), see figure 
1191 The spike includes a close pair with a separation large 
enough to be individually spatially resolved in the spectrum. 
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Figure 16. Field J1040. 7-1155: The redshift distribution of grade 
A confirmed galaxies is shown by the dashed red hne and the 
distribution of all confirmed galaxies (grade A and B) is shown 
by the dotted blue line binned with the same parameters as the 
individual field data. The solid green line is the total redshift 
distribution averaged over the ten survey fields. The sky emission 
spectrum offset from zero at the redshift of Lyo is shown by the 
solid black line. 
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Figure 18. One dimensional spectra of two LBGs separated by 
1.1" or 7kpc on the sky. The sky spectrum at this redshift is 
plotted below each spectrum and the location of Lya at z = 4.951 
is shown. 



This pair was identified as a Lya emitter and break-only 
galaxy as shown in figure 1181 When both grades of spectro- 
scopic confirmations of high redshift galaxies are included, 
this overdensity contains nine galaxies. If we include objects 
between z — 4.95 and z = 5.15, i.e. Sz = 0.2, 17 grade 
A and B galaxies are included. However this redshift range 
is much larger than that expected for a cluster of galax- 
ies or any other bound structure. The broad nature of this 
spike may suggest we are observing a filamentary structure 
end-on. Again the overdensity is mainly determined by line 
emitters, a reflection of the greater fraction of line emitting 
galaxies confirmed compared to break-only galaxies in this 
field. 

The galaxies in the redshift spike are spatially dis- 
tributed over a smaller fraction of the field than those in 
J1040. 7-1155. They are clustered to one side of the 7 x 7 
arcminute field, as shown in figure [20] The closest grade A 
galaxy pair has a co-moving separation of 670 kpc (excluding 
the very close pair), with the overdensity of galaxies cover- 
ing 27Mpc (co- moving distance between furthest grade A 
confirmed galaxies in three dimensions). Again, the spectra, 
colours and morphologies of the galaxies within the redshift 
spike are typical of the sample as a whole. 



Figure 17. Spatial distribution of spectroscopically confirmed 
grade A galaxies (diamonds) and grade B galaxies (crosses), which 
populate the redshift spike, across the field J1040. 7-1155. Coor- 
dinates are rel ative to the fie l d centr e. The asterisk denotes the 
position of the lStanwav et al.l l l2008al ) CO source discussed in the 
text. 



8.2 The Nature of the Large Scale Structure 

As discussed above, two of the ten survey fields show evi- 
dence of large scale structure in the form of spikes in their 
redshift distributions. Such overdensities could be a forming 
cluster, sheet or filament being probed in these fields. Given 
that across the survey as a whole zero to one sources per 
0.1 in redshift are expected, the identification of two fields 
with very clear overdensities of many times this expected 
number is highly significant. However, given the stochastic 
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Figure 19. Field J1054. 7-1245: The rcdsiiift distribution of grade 
A confirmed galaxies is shown by the dashed red line and the 
distribution of all confirmed galaxies (grade A and B) is shown 
by the dotted blue line. The solid green line is the total redshift 
distribution averaged over the ten survey fields binned with the 
same parameters as the individual field data. The sky emission 
spectrum offset from zero at the redshift of Lyo is shown by the 
solid black line. 
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Figure 20. Spatial distribution of spectroscopically confirmed 
grade A galaxies (diamonds) and grade B galaxies (crosses), which 
populate the two 0.1 redshift spikes (distinguished by bold sym- 
bols), across the field J1054-7.1245. Coordinates are relative to 
the field centre. 



properties of LBGs (|Verma et al.ll2007f ) , the nature of these 
overdensities is not immediately clear. 

As the typical star-formation lifetime of high redshift 
LBGs is a few tens Myrs (assuming the ma jority of the 
galax ies within the overdensities are young, IVerma et al.l 
[20071) , the apparent coordination or synchronisation of these 
short-lived star-bursts over the entire field needs explana- 
tion. Given the distances between the galaxies within the 
redshift overdensities, it is highly unlikely that these galaxies 
are dynamically bound or interacting sufficiently to have a 
common starformation trig ger. Assuming the y oung ages of 
these galaxies predicted by IVerma et al.l (|2007|), <45Myrs, 
it is likely that the galaxies have only travelled 50kpc, out- 
side the Hubble flow, in their lifetime (assuming the galaxy 
is travelling at 1000 kms"'^). This suggests that if two galax- 
ies are more than 100 kpc apart, it is unlikely that they 
have interacted. With the closest pair of galaxies (in fleld 
J1054. 7-1245) separated by 670 kpc (co-moving distance in 
three dimensions), the star formation observed cannot be 
connected. This limited interaction region is reflected in the 
high resolution imaging of the HST data. It was found that 
many objects detected in the ground based data were in fact 
multiple systems, at most 16 kpc apart and are subsequently 
close enough to potentially have star formation triggered by 
the mutual interaction. 

The apparent synchronisation of the relatively short- 
lived starbursts is likely to be an illusion for the above rea- 
sons. A more plausible explanation for the overdensity of 
objects undergoing a relatively short-lived phenomenon is 
that there is far more baryonic material in the large-scale 
structure that hosts these galaxies than is visible in the rest- 
frame UV. At any one time a small fraction of this material 
is visible as unobscured UV-bright starburst episodes lasting 
a few tens of Myrs. If one were to wait a few tens of millions 
of years and re-observe the region, the large scale structure 
would still be marked out by an overdensity of UV-luminous 
starbursts, but probably different ones to those seen today. 
Therefore, there is no synchronisation of the starbursts, just 
a high probability that several regions of the large-scale 
structure will be 'on' at any one time. One way to try to 
detect this 'dark' material connecting the observed systems, 
is to look for line or continuum emission, such as carbon or 
CO emission lines, from gas within the structure, in the mil- 
limetre and submillimetre regimes. Such an experiment was 
carried out in a 2 arcminute diameter region centred between 
two L BGs with concurrent redshifts in one of the overdense 
fields (|Stanwav et al.ll2008al ) . Although none of the optically 
identified LBGs were detected in the 37.7 GHz band, a faint 
C0(2-l) emission line was detected at the same redshift as 
the optically identified overdensity (and subsequently recon- 
firmed in independent data, Stanway et al., in prep), shown 
in figure [T71 With no optical or near-infrared counterpart in 
the available imaging and a inferred molecular hydrogen gas 
mass of 2 X 10^" M0 this line emitter is tracing previously 
undetectable baryonic mass. 

The assumption that there is a large reservoir of un- 
detected nM±erial_Ja_Hie_structures is consistent with the 
work of iThomas et al.l (|2005|), which, through analysis of 
fossil populations in nearby galaxes, determined that mas- 
sive galaxies found in high density environments formed 
the majority of their stars at 2 > 4. Further support for 
such old star formation in massive galaxies is provided by 
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the distant r ed ga laxy (DRG) population at 2 = 2 — 3. 
iLabbe et all (|2005h observed DRGs which had UV to 8 ^m 
spectral energy distributions fitted best by an old stellar 
pop ulation with very l ittle c urrent star formation. From 
the lBruzual fc Charloj l|2003l ) models used, they infer that 
the majority of the stellar mass was in place by z ~ 5 or 
earlier (depending on the metallicity assumed). Similarly, 
IPanter et al.l ((20071) report that the most massive systems 
seen today formed the majority of their stars at z > 4, 
from analysis o f 30,000 galaxy sp ectra from the Sloan Digi- 
tal Sky Survey. IStark et al.l (|2007l ) also predict that a signif- 
icant fraction of the mass at 2; ~ 5 is not being observed in 
Lyman break galaxy surveys. Assuming that galaxies seen 
here at z ^ 5 are m a inly y oung, the typical stellar mass 
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found bv lVerma et al.1 (|2007l l was a few 10^ M©. If this mass 
is applied to the galaxies in the redshift overdensities, and 
it is assumed that over time they will merge into one of the 
most massive systems today (a best-case scenario to accu- 
mulate the most mass), the observed stellar mass at z = 5 
would only account for around 2% of the final baryonic mass 
of a massive galaxy at z = 0. Although the structures may 
extend beyond the 44 arcmin'^ field of view of the imaging 
or a fraction of the galaxies could be older and more massive 
than typical of the population, there is still large fraction of 
the mass unaccounted for. 

The apparent problem of starburst synchronisation can 
be avoided if the galaxies in the overdensities have older stel- 
lar populations with longer-lived starbursts than are typical 
of the population as a whole. Alth o ugh t he majority of LBGs 
at z ~ 5 studied by IVerma et al.l (|2007l ) are young, a num- 
ber of older, more massive systems were identified. Similarly, 
individua l relatively brigh t obj ects at z ~ 6 have been ob- 
served bv lYan et all (|2006l ) and lEvles et al.l (|2005l 'l and were 
found to have older and more massive stellar populations. 
Finding individual old galaxies suggests that the formation 
and assembly of massive galaxies had already started by 
z ~ 10. As such sources are expected to be in the minority 
of the z ~ 5 population, discovering a connected system of 
old galaxies would pose significant problems for our under- 
standing of galaxy and structure formation. 

To distinguish between these two scenarios, multiband 
SED fitting of the sources is required. With the additional 
deeper near infrared and Spitzer IRAC imaging that is cur- 
rently being added to this extensive data set, it will be pos- 
sible to probe the rest-frame UV to optical spectral energy 
distributions of these galaxies and determine whether they 
have young or old stellar populations. 



9 Z - 6 GALAXIES 

As noted previously in section [2]2j the imaging dataset of V, 
R, I, z, J and A's-bands can be used to select high redshift 
objects at z; ~ 6 using an Iab — zab colour cut. A total 
of 73 candidates were observed spectroscopically with an 
Iab — Zab > 1-3 and non detections in the V and _R-bands. 
From this sample two objects were confirmed to be galaxies 
at z ^ 6. 

A galaxy with an observed zab ~ 25.7 and a colour of 
Iab — Zab = 1.4 was confirmed to be at z = 5.99 with strong 
line emission and weak continuum at longer wavelengths but 
none at shorter wavelengths. With a Lya fiux of 1.9 x 10~^^ 



Figure 21. 2D spectra of the LBGs a,t z = 5.99 (top frame) and 
z = 6.41 (bottom frame). Both show the spectral region between 
7800A and 9800A. 



ergs cm"'^ s~^, the rest-frame equivalent width is 46A mea- 
sured using the same method as the z ~ 5 galaxies (section 
[S| . The second galaxy was confirmed with a weaker Lya line 
and continuum break at z — 6.41. This galaxy was observed 
to have zab = 24.9 and a colour of Iab — zab > 2. The 
line fiux and rest-frame equivalent width of this emission 
are 9.4 x 10~^* ergs cm"'^ s~^ and 8A respectively. 

Of the other candidates observed, 8 were at low redshift 
and the remaining 63 had signal-to-noise too low to be able 
to identify the sources. The large number of unidentified 
objects is unsurprising given the decreasing response of the 
FORS2 instrument beyond 9000A. With better instrumen- 
tal responses and a wider range of wavelengths uncluttered 
by sky lines at shorter wavelengths, it is markedly easier 
to study z ~ 5 galaxies than their higher redshift coun- 
terparts, even those observed when the Universe was just 
200-250 Myrs older. 



10 CONCLUSIONS 

We have carried out a large spectroscopic survey of 10 
widely-separate fields targeting z --^ 5 LBGs. The candi- 
dates were selected using deep optical imaging combined 
with a colour selection tuned to identify the Lyman Break 
at z: ~ 5. A total of 70 LBGs at 4.6 < z < 5.6 have been 
identified, 38 galaxies with Lya emission ranging in strength 
from 2.6 x 10~^* to 7 x 10"^''' erg cm~^ s~^ (luminosities 
between 10*^ and 10'*'^ erg s~^) , and 32 were identified by 
their strong continuum break at a rest wavelength of 1216A. 
An additional two galaxies were confirmed to be at z; ~ 6, se- 
lected for their red l — z colours. A comparison of the equiva- 
lent width distribution of the Lya emission of the confirmed 
z_~5 LBGs and that of z ~ 3 LBGs from IShaplev et al.l 
([2003) shows that the confirmed higher redshift objects are 
twice as likely to have large (> 20 A) equivalent width Lya 
emission. Just over half of the z ~ 5 objects have Lya lines 
stronger than this. The largest measured equivalent width is 
Wo = 122±18A. However, including plausible corrections for 
the spectroscopically-unconfirmed objects which are never- 
theless at these redshifts removes any significant difference in 
the fraction of sources with large (> 20A) equivalent widths 
in both samples. 

A comparison of the average optical spectra and av- 
erage Iab — Kab colours of LBGs with and without de- 
tectable Lya emission showed that the break-only galax- 
ies had a redder colour and hence steeper rest-frame UV 
continuum slope than the line emitting galaxies. It is dif- 
ficult for stellar age differences alone to explain this dif- 
ference given the statistics on the ages of sources from 
IVerma et al.l (|2007l ). It can be straightforwardly explained 
by variati ons in reddening (a d ifference of E(B-V)'^ 0.2 as- 
suming a ICalzetti et al.l (|200Gr ) reddening law can account 
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for the difference in colours) and/or metallicity between the 
two subsamples, though age could clearly also play a part. A 
similar, but smaller colour difference is seen between samples 
of line-emitting and break-only LBGs at lower redshifts. 

If the local correlation between UV-slope and metallic- 
ity is assumed to be applicable at high redshifts, a typical 
z ^ 5 LBG has a metallicity of approximately 0.1-0.2 of the 
solar value, approximately a third that of the metallicity of 
2 ~ 3 LBGs. 

The morphology of the confirmed z ~ 5 galaxies was 
studied using high resolution HST/ACS images and was 
found to be similar to that seen in larger photometric sam- 
ples with a typical half-light radius of 0.15" or Ikpc. How- 
ever, a higher fraction of multiple and disturbed systems 
was observed (16 out of 18 LBGs) suggesting such systems 
are more likely to be confirmed spectroscopically, with a 
greater incidence of Lya emission. A comparison of the mul- 
tiple z ~ 5 systems with local analogues of LBGs suggests 
that the components identified in the high resolution imag- 
ing could be either super-starburst regions within a larger 
structure with a surface brightness below the detection limit, 
or part of merging objects. 

The surface densities of LBGs in two of the ten fields 
are substantially higher than in the rest. In both cases the 
fields contain a clear excess of LBGs over a narrow range 
in redshift {z = 5.05 ± 0.1 and z = 5.15 ± 0.05). Although 
these are coherent structures, they may extend beyond the 
observed fields (with scale sizes of 20-50 Mpc) and cannot 
be bound or virialised systems. Given the distances between 
the LBGs within these excesses and the expected short (no 
more than a few tens of Myr) timescale for typical z ^ ^ 
LBGs, the excess of systems cannot be caused by gravita- 
tional interactions between the LBGs triggering starbursts. 
It is more likely that the LBGs are comparatively rare UV- 
luminous markers for a much larger mass (both baryonic 
and dark) of UV-dark material in these volumes. 

In addition to the survey for z ~ 5 LBGs, we use the 
same photometric data to carry out a comparatively limited 
search for z ~ 6 LBGs by looking for objects with spectral 
breaks between the / and z-band images. We confirmed two 
such galaxies, one at z = 5.99 with a Lya line flux of 1.9 x 



10" 



erg cm '^ s ^ and zab ~ 25.7 and another at z — 6.41 



^* erg cm ■^ s 



with a much weaker Lya line (9.4 x 10 

but considerably brighter continuum {zab = 24.9) 
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